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We show that the dominant Gamow-TeUer part, AfgTi oi the nuclear matrix element governing 
the neutrinoless /3/3 decay is related to the matrix element M^l' governing the allowed two-neutrino 
/3/3 decay. That relation is revealed when these matrix elements are expressed as functions of the 
relative distance r between the pair of neutrons that are transformed into a pair of protons in the 
/3/3 decay. Analyzing this relation allows us to understand the contrasting behavior of these matrix 
elements when A and Z is changed; while Mq'j, changes slowly and smoothly, M^" has pronounced 
shell effects. We also discuss the possibility of phenomenological determination of the M^l' and from 
them of the Mq^ values from the experimental study of the strength functions. 



I. INTRODUCTION 

Observing 0i/f3/3 decay would tell us that the to- 
tal lepton number is not a conserved quantity, and 
that neutrinos are massive Majorana fermions. An- 
swering these questions is obviously a crucial part 
of the search for the "Physics Beyond the Standard 
Model" . Consequently, experimental searches for 
the QvPj3 decay, of ever increasing sensitivity, are 
pursued worldwide (for a recent review of the field, 
see e.g. ^). However, interpreting existing results 
as a determination of the neutrino effective mass, 
and planning new experiments, is impossible with- 
out the knowledge of the corresponding nuclear ma- 
trix elements. Their determination, and a realistic 
estimate of their uncertainty, are therefore an inte- 
gral part of the problem. 

The nuclear matrix elements M^'^ of the Oi^/3/3 de- 
cay must be evaluated using tools of nuclear struc- 
ture theory. Unfortunately, there are no observ- 
ables that could be simply and directly linked to 
the magnitude of Oi^/3/3 nuclear matrix elements 
and that could be used to determine them in an 
essentially model independent way. In the past. 
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knowledge of the 2z//3/3-decay rate, and therefore 
of the corresponding matrix elements M^*^, and of 
the ordinary /3 decay ft values and the correspond- 
ing beta strength distributions, were used to con- 
strain the nuclear model parameters, in particular 
when the Quasiparticle Random Phase Approxima- 
tion (QRPA) was employed 12H5]. In the present 
paper we discuss a novel relation between these nu- 
clear matrix elements. 

Very early, Primakoff and Rosen I6j speculated 
that since the operators governing Oi^/3/3 and 2vf3[} 
decays differ by a relatively gentle radial depen- 
dence, approximately of the form 1/rij, the corre- 
sponding matrix elements might be proportional to 
each other with the proportionality constant ^ 1/R, 
where R = 1.2^^/^fm is the nuclear radius. At that 
time the authors also believed that the 21/(3/3 de- 
cay can be treated in closure, thus that the corre- 
sponding matrix element is dimensionless, while in 
fact the realistic 2i>/3/3 matrix element has dimen- 
sion energy"^. Also, the 01^/3(3 matrix elements are 
now, by convention, made dimensionless by includ- 
ing the nuclear radius i? as a multiplicative factor, 
which is compensated by the factor R^^ in the cor- 
responding phase space function G^"^. 

Modern nuclear structure evaluations of these ma- 
trix elements do not support the conjecture of pro- 
portionality between M"'^ and M'^" . The rate of the 
2i//3f3 decay has been determined experimentally in 
many nuclei, and hence the 2v(3P matrix elements 
M'^'^ are known. They exhibit pronounced shell ef- 
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fects and vary rather abruptly between nuclei with 
different Z and A. At the same time, the calcu- 
lated Ov/if] nuclear matrix elements, whether based 
on the QRPA |S{5], nuclear shell model [7HS], or the 
Interacting Boson Model |T0], do not show such a 
variability; instead they vary relatively smoothly be- 
tween nuclei with different Z, A. The reason for the 
difference is, presumably, the very different momen- 
tum transfer q involved in these matrix elements, 
even though they involve the same initial and final 
nuclear states. In the 2iy/3/3 decay the momentum 
transfer q is restricted to g < Q, where Q is the en- 
ergy difference of the initial and final atomic masses. 
Hence, the allowed approximation is valid, qR <C 1, 
and only the Gamow- Teller operator crT+ and only 
the 1+ virtual intermediate states, contribute. On 
the other hand, in the Qvj3j3 decay the momentum 
transfer is of the order of the nucleon Fermi momen- 
tum q ^ 200 MeV, qR> I, and all J'" virtual in- 
termediate states can contribute significantly. Our 
discussion here sheds more light on the different be- 
havior of the Oi^/3/3 and 2i//3/3 matrix elements. 

It is worthwhile to remember another type of re- 
lation, explored in the classic paper by Pontecorvo 
[TT] . At that time the available information on /3f3 
decay was based on the geochemical determination 
of the total decay rate, l/rtot = ^/to^ + ^/t2i^- Since 
these two modes scale very differently with Q (~ 
for and ^ Q^^ for 2i/) Pontecorvo suggested that 
comparing the total lifetimes of two isotopes, ^•^°Te 
and ^^^Te, which have very different Q values, might 
reveal the presence of the lepton number violating 
Ov decay, provided the nuclear matrix elements of 
these two isotopes are identical. While the matrix 
elements of these two isotopes are indeed rather 
close, they are not quite the same. Moreover, we 
know today that the Ov decay rate is very much 
smaller, if it is indeed nonvanishing, than the 2^ 
decay rate. 

The present paper is structured as follows. In the 
next section II we describe the formalism that leads 
to the relation between the Gamow- Teller part of 
the Qvj3l3 matrix element and the 2i//3/3 matrix el- 
ement evaluated in the closure approximation. We 
also discuss the validity of the closure approxima- 
tion in the Oi^/3/3 case. In the following section III 
we discuss this novel relation in more detail and 
show numerous examples. In section IV we briefly 
discuss the issue of quenching of the axial current 
matrix elements. While closure is a rather poor ap- 
proximation in the 2i//3/3 case, we argue in section V 
that combining the known lifetimes with the often 
measured distribution of the P~ and strengths 
constrains the M'^f values substantially. We believe 
that the relation found here allows one to better 
understand the different behavior of these matrix 
elements. We conclude in the last section. 



II. FORMALISM 

Assuming that the Qu(5(i decay is caused by 
the exchange of the light Majorana neutrinos, the 
halflife and the nuclear matrix element are related 
through 
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Tl/2 



^G'-^{Q,Z)\M''^\' \{fnpp)\' , (1) 



where G^'^{Q, Z) is the easily calculable phase space 
factor, {mpp) is the effective neutrino Majorana 
mass whose determination is the ultimate goal of 
the experiments, and M^'^ is the nuclear matrix ele- 
ment consisting of Gamow- Teller, Fermi and Tensor 
parts. 



Ml" 



AfO- = Mg^j,(l + XF + XT) 



(2) 

where xf and xt are the matrix element ratios that 
are smaller than unity and, presumably, less depen- 
dent on the details of the applied nuclear model. In 
the following we concetrate on the GT part, M^.^, 
which can be somewhat symbolically written as 



(3) 



where -ff(r/fc,i?) is the neutrino potential described 
in detail below and rj^ is the relative distance be- 
tween the two neutrons that are trasformed in the 
decay into the two protons. 

In Ref. |4], based on the QRPA, as well as in Ref. 
[7] based on the nuclear shell model, the function 
C^"(r) that describes the dependence of the M^'" 
on the distance r/^ was introduced. Formally, this 
function can be defined as [12] 



(4) 

where ^(x) is the Dirac delta function. Obviously, 
this function is normalized by 



{r)dr , 



(5) 



and has the dimension lenght ^. The shape of 
C^xir) is very similar in both QRPA and NSM and 
in all cases consists of a peak with maximum at 
r ^ 1 fm ending near r ^ 2.5 fm, and of very little 
contributions for larger values of r. 

Now lets turn to the case of the 2z/ decay mode. 
The matrix element M^" governing the 2i//3f3 decay 
mode is of the form 



+ ||m)(TO||(Tr"'"||i) 
(Mi + Mf)/2 



(6) 
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FIG. 1: Multipole decomposition of C^i{r) as function 
of relative distance of two /3-decaying neutrons in the 
''^Ge nucleus. Calculation performed for ^^Ge with 23 
single particle levels model space. Positive parity mul- 
tipoles are shown in the upper panel and the negative 
parity ones in the lower panel, (color online) 

where the sumation extends over all 1+ virtual in- 
termediate states. We can introduce also the clo- 
sure analog of A'P" , denoted by M'^i , by replacing 
the energies Em by a properly defined average value 
E^v Thus, 

M^l' = Af2. ^ (^^^ _ ^ ^,^^)/2) . (7) 

In analogy with Eq. Q we can define the new func- 
tion 

poo 

Mir = / Clr{r)dr . (8) 
Jo 

While the matrix elements A/^'' and M^f get con- 
tributions only from the l"*" intermediate states, the 
function C^i gets contributions from all interme- 
diate multipoles. This is the consequence of the 



FIG. 2: C^i'{r) as a function of the relative distance 
of the decaying neutron pair for different nuclei, (color 
online) 

5 function in the definition of Cli{r). When ex- 
panded, all multipoles contribute. Naturally, when 
integrated over r only the contributions from the 
1"*" are nonvanishing. An example of the multipole 
decomposition of Cli{r) is shown in Fig. [ij and in 
Fig. plwe show the functions Cliijr) for a variety of 
/3/3 decaying nuclei. 

For completeness we show here the QRPA formula 
used for the evaluation of the function Cli{r) and its 
multipole decomposition depicted in Fig. [l] First, 
the function 

/„^„',P.,'M= (9) 
{p{l),p'{2){r)-J\\ a^-cj2 II n{l),n'{2){r)-J) 

is introduced where r is the relative distance be- 
tween the neutrons in the states n and n', respec- 
tively protons in p and p' . Then, the part of Cli{r) 
with the multipolarity J'^ is given by 

-Mj^J pnp'n' 
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xy2jTT X I J I X /„^„, p,^,(r) X (10) 

{o+\\[cpnMJ''kf){rkf\rh){rkM4cn]j\\ot). 

Here fc^ and kf are the labels of the excited 
states with the niultipolarity J'^ in the interme- 
diate nucleus built on the initial and final nu- 
clear ground states, and {0^\\[Cp,Cn']j\\J'^kf) and 

( J'^ fci 1 1 [Cp c„],/ 1 10^) are the corresponding QRPA 
amplitudes. 

It is now clear that, by construction, 

C°Jr{r)^H{r,E)xCf{r), (11) 

which is valid for any shape of the neutrino poten- 
tial H{r, E). Thus, if C^l'lr) is known, Cg'^(r) and 
therefore also Mq'j, can be easily determined. The 
equation (111 represents the basic relation between 



the Ov and /3/3-decay modes that we will explore 
further. 

Note that while the function C^/'(r) has a sub- 
stantial negative tail past r ~ 2 — 3 fm, these dis- 
tances contribute very little to CQ^(r). This is a 
consequence of the shape of the neutrino potential 
H{r, E) that decreases fast with increasing values of 
the distance r. 



A. Neutrino potential 

The neutrino potential HQT{r, E) governing the 
Gamow- Teller part of the matrix element M^'^ is 
defined as 



2R 

■K 

where 



i?GT(r,So.) = 

jo{qr) ^ fFNsi<l^)9HOT{q^)dq ,(12) 
q + Eoi, 



fpNS 
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(13) 



takes into account the finite size of the nucleon and 
is usually approximated using the above dipole type 
form factor with Ma = 1.09 GeV [T^ (varying Ma 
between 1.0-1.2 GeV makes little difference). The 
function gHOxiq^) includes the terms from higher 
order hadron currents, namely induced pseudoscalar 
and weak-magnetism |14| . The short range correla- 
tions are included using the method of Ref. [T5] . 
The Jastrow-like two-body function derived there is 
applied when the radial integrals in both functions 
C"'^ and C^l' are evaluated; they do not appear ex- 
plicitly in eq. (12). 



We show in Fig. [3] the shape of the potential. 
When the finite nucleon size, higher order terms are 



neglected, and Eq^ = is assumed, the potential 
has Coulomb-like shape R/r. The full potential, Eq. 
(12 1, however, is finite at r = 0, H{r — )■ 0,Eq^ = 
0)"= 5MaR/IG. Including the higher order currents 
and finite E in Eq. ( 12 1 increases the value of H{r — 
0) by -30%. 
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FIG. 3: The potential _ffGT(?', -B). Diflterent approximate 
forms, as well as the exact one, are shown (color online). 
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FIG. 4: Matrix elements Af"" for the indicated nuclei 
evaluated in the closure approximation as a function of 
the assumed average excitation energy (color onhne). 
The values of A/"" obtained without the closure approx- 
imation are 5.24C^''Ge), 2.62 C'^^Zr), 4.99 (i°°Mo), and 
4.07 (^^°Te). 



B. Validity of the closure approximation for 
the OvPP matrix element 

The closure approximation, i.e. the replacement 
of the summation over the virtual intermediate 
states by matrix element of a two-body operator, 
is typically used in the evaluation of the M'^" . It 
is worthwhile to test the validity of this approxi- 
mation. Such test can be conveniently performed 
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within the QRPA, where the sum over the inter- 
mediate states can be easily carried out. In fact, 
the calculations performed in Refs. [2H1] do not use 
closure. In this context one can ask two questions: 
How good is the closure approximation? And what 
is the value of the corresponding average energy? 
In Fig. |4] we illustrate the answers to these ques- 
tions. The exact QRPA matrix elements shown in 
the caption can be compared with the curves ob- 
tained by replacing all intermediate energies with a 
constant E, which is varied there between and 12 
MeV. One can see, first of all, that the M'^^ changes 
modestly, by less than 10% when E is varied and, 
at the same time, that the exact results are quite 
close, but somewhat larger, than the closure ones. 
Thus, using the closure approximation is appropri- 
ate for the evaluation of M'^'^ even though it slightly 
underestimates the M"*^ values. However, the cor- 
responding uncertainty is not more than the other 
uncertainties involved. We compare the QRPA ex- 
act and closure M^'^ for all nuclei of interest in the 
next section. 



decay half-lives of Ref. [17] . 

In Table |l] we show both the calculated Oi^PfS 
NMEs evaluated with and without the closure ap- 
proximation, as well as only the GT parts of their 
values. Also shown are the experimental 2v/3f3- 
decay NMEs. Using the QRPA method the clo- 
sure matrix elements M^l" were also evaluated. One 
can see that the spread among the candidate nuclei 
of the 21/(3(3 NMEs is significatly larger when com- 
pared with the spread of the calculated 0iy(3(3-decay 
NMEs. The table also demonstrates that using the 
closure approximation for evaluation of A/"'' makes 
relatively little difference and that the GT part of 
M°'' is dominant in all considered nuclei. 

The values of M^" in Table | might be com- 
pared with the corresponding entries in Table II of 
Ref. [iS] . There are small differences between them 
caused by several changes made in the present work. 
We use now the updated values of T^^ of Ref. [T7] 
and the more realistic gA = 1.269 instead of 1.25. 
In evaluating M^"^ we adjust here the energy de- 
nominators such that the first 1+ state has the ex- 
perimentally known energy value. Moreover, the 
present results are based on the level scheme with 
just 23 single particle states, while Ref. [T5] uses an 
average of several sets of single particle energies. 

Another characteristic feature of the relation be- 
tween the M^^rp and M^l" is illustrated in Fig. [H] 
There we show the integrals, i.e. the functions of 
the upper limit of the integration. 
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III. RESULTS AND DISCUSSION 

We evaluated the nuclear matrix elements (NME) 
Mo^rp and M^'^ with and without closure approxi- 
mation using the Quasiparticle Random Phase Ap- 
proximation (QRPA). 

For all nuclear systems the single-particle model 
space consisted of O—dfiLu oscillator shells plus Oiii/2 
and 0ii3/2 levels both for protons and neutrons (23 
single particle states). The single particle energies 
were obtained from the Coulomb-corrected Woods- 
Saxon potential. Two-body interaction G-matrix 
elements were derived from the Argonne V18 one- 
boson exchange potential within the Brueckner the- 
ory. The pairing interaction was adjusted to fit the 
empirical pairing gaps |16j . The particle-particle 
and particle-hole channels of the G-matrix inter- 
action of the nuclear Hamiltonian H were renor- 
malized by introducing the parameters gpp and gph, 
respectively. While gph = 1.0 was used through- 
out, the particle-particle strength parameter gpp 
was fixed by the data on the two-neutrino double 
beta decay rates [5H1] for each nucleus separately. 
In the calculation of the Oi^/3/3-decay NMEs the 
two-nucleon short-range correlations derived from 
same potential as residual interactions, namely from 
the Argonne potential |15) . were applied. The un- 
quenched value of the axial current coupling con- 
stant, gA = 1.269 was used here. The modifications 
caused by the quenching of the weak axial current 
are discussed in the following two sections. 

On the other hand, the absolute values of M^^^ 
were deduced from the averaged values of 2v(3(3- [H] demonstrate that the matrix elements Mq^j. are 



Obviously, 

As one can see the integrals lorifo) saturate for 
To > 2 — 3 fm since the function CQ^(r) is very 
small past these values of r. On the other hand, the 
functions l'^'^ {tq) change drastically, even becoming 
sometimes negative, for > 2 — 3 fm. That is a 
reflection of the behavior of the function C'^i (r) that 
has a substantial tail for rp > 2 — 3 fm. In addition. 
Fig. [5] also demonstrates that the corresponding 
integrals are almost independent on the number of 
included single-particle states, as long as at least 
two full oscillator shells are taken into account. 

Remembering that in a nucleus the average dis- 
tance between nucleons is ^1.2 fm we can somewhat 
schematically separate the range of the variable r 
in the functions C^'^(r) and Cli{r) into the region 
r < 2-3 fm governed by the nucleon-nucleon corre- 
lations, while the region r > 2-3 fm is governed by 
nuclear many-body physics. The integrals in Fig, 
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TABLE L The nuclear matrix elements IMq^I and IM""] calculated within the QRPA. For the param- 

eters used, see the text. The 2i'/3/3-decay nuclear matrix element |Mg^p| were deduced from the avaraged values of 
the 2!//3/3-decay half-lifes [H] and Mf," were obtained within the QRPA. For ^^^Xe, where only the upper limit of the 
2v half-life exists, the range shown covers the range of half-lives from the experimental limit to infinity. All entries 
are evaluated with qa = 1.269. 

NME '""Ge ^^ge ""^Zr ^°°Mo ^^''Cd ^^^Te '^°Te '^"Xe 

2i//3/3-decay NMEs 
0.136 0.095 0.090 0.231 0.126 0.046 0.033 (0,0.031) 
0.099 -0.126 -0.802 -0.933 0.059 -0.462 -0.464 (-0.41 -0.25) 

0!//3/3-decay NMEs within closure approximation 
4.12 3.61 1.89 3.72 2.77 3.63 3.09 (1.61,1.83) 
5.02 4.44 2.34 4.59 3.36 4.44 3.79 (2.00,2.24) 

Oi/PjS-decay NMEs without closure approximation 
4.33 3.82 2.16 4.10 2.91 3.92 3.36 (1.76,1.96) 
5.24 4.65 2.61 4.99 3.51 4.75 4.07 (2.15,2.38) 
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FIG. 5: Integrals /^"(ro) and /St(7-o), Eq.(|l4|, as func- 
tion of the upper limit ro. Three different spaces of 
single-particle states are considered, small, medium and 
large (color online). 



almost independent of the "nuclear" region of r and 
hence one does not expect rapid variations of their 
value when A or Z of the nucleus is changed. On 
the other hand, the 2^ closure matrix elements de- 
pend sensitively on that region of r and hence one 



expects sizable shell effects, i.e. a significant varia- 
tion of M^'^ and M'^l' with A and Z, in agreement 
with observations. 



IV. QUENCHING OF THE AXIAL 
CURRENT MATRIX ELEMENTS 

It is well known that Gamow- Teller /3-decay 
transitions to individual final states are notice- 
ably weaker than the theory predicts. That phe- 
nomenon is known as the axial current matrix ele- 
ments quenching. The /3-strength functions can be 
studied also with the charge exchange nuclear reac- 
tions and a similar effect is observed as well. Thus, 
in order to describe the matrix elements of the oper- 
ator err, the empirical rule (crT)g^j ~ 0.6(crT)^g^g; 
is usually used (see [T5H2D]). Since these operators 
accompany weak axial current, it is convenient to 
account for such quenching by using an effective cou- 
pling constant g"^^ ^1.0 instead of the true value 
gA = 1.269. 

The evidence for quenching is restricted so far to 
the Gamow- Teller operator ctt and relatively low- 
lying final states. It is not known whether the other 
multipole operators associated with the weak axial 
current should be quenched as well. In fact, the 
analysis of the muon capture rates in Refs.[lTJ [52] 
suggests that quenching is not needed for this pro- 
cess with momentum transfer q ^ 100 MeV. 

Since the 21^/3/3 decay involves only the GT opera- 
tors and relatively low-lying intermediate states, one 
could expect that the quenching might be involved 
in that case. Whether it should be included also 
for the Ov mode remains an open question. In the 
previous paper, Ref. [3], it was shown that by mak- 



7 



ing the adjustment of the particle-particle coupling 
strength gpp so that the experimental 2^ halflives are 
correctly reproduced, the predicted Ov decay rates 
are affected by the possible quenching less than the 
ratio [{g'^'^ /gA]'^ might suggest. 

Following Ref. [3] we define the "quenched" nu- 
clear matrix elements 



1.269 



M^'^ig', 



and use the analogous definitions for M'J'^ 
and for the integral /'^"^(ro) 
p). 



and 



I'Sriro) 



(15) 



see eq. 



We use this definition since the experimental 
quantities, the halflives T^'^ and T^'', are then sim- 
ply proportional to 1/M'^ without the need to mod- 
ify the phase space factors G^^ or G^"^ when a dif- 
ferent value of g'^-^ is used. Note that as a con- 
sequence of our choice of renormalization of the 
particle-particle coupling constant gpp the 2vl3j3 ma- 
trix elements M"^^ by definition become independent 
of g"/^ and thus M"^" = M'^"^. 

In Fig. ^ we show the integrals I'^'^{ro) and 
^Gt('"o) for the case of the decay of ™Ge and three 
values of {g'^^)- One can see that in the case of 
M'J'^ not only does the final value depend on {g'^/^), 
but it affects the dependence on the distance tq as 
well. With the standard gA = 1.27 the peak at 
tq ^ 2 fni is compensated by the long tail extending 
to much larger tq, while for the heavily quenched 
s'a^ = 0.8 the I'^^iro) almost saturates at the 
much smaller values of rp. In contrast, the integrals 
^Gt('^o) saturate at tq ^ 3 fm for all considered 
values of g'^'^ ■ 

It was shown in Ref. [3] that with g'^-^ — 1.0 the 
full matrix elements Al"^'^ are reduced by 10-15% 
compared to their value with gA = 1.25 used in that 
paper. Here we use the more correct gA = 1.269, 
for adjustment of the particle-particle coupling con- 
stant gpp we use the 2v halflives of Ref. that 
differ slightly from the halflives used in [3], and 
the treatment of the short-range correlations is now 
based on the Ref. [15 while in 3 it was based 
on the phenomenological Jastrow-type function of 
Ref. [23]. In the present work the matrix elements 
M'^" are 20-30% smaller with g^/^ = 1.0 than with 
gA — 1.269. Similar effects are also visible in Table 
II of Ref. HSl. 



DETERMINATION OF THE MATRIX 
ELEMENT Mil" 



While the nuclear matrix elements M are sim- 
ply related to the 2v half- life Ty2^ and are therefore 
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FIG. 6: The running sums of I'^^ (upper panel) and 
^GT (lower panel) for ^^Ge and different effective values 
of qa- (color online). 



known for the nuclei in which Tfj^ has been mea- 
sured, the closure matrix elements need be de- 
termined separately. There are several ways how to 
accomplish this task: 

1. Rely on a nuclear model (e.g. QRPA or nu- 
clear shell model), adjust parameters in such a 
way that the experimental value of M"^^ is cor- 
rectly reproduced, and use the model to eval- 
uate M'^i . (In QRPA the usual adjustment 
is the renormalization of the particle-particle 
coupling constant gpp so that the T'^j2 '-'-'r- 
rectly reproduced.) This procedure is used in 
Table H 

2. Use the measured /3~ and strength func- 
tions and assume coherence (i.e. same signs) 
among states with noticeable strengths in 
both channels. In this way an upper limit of 
M^j" can be obtained. 

3. Finally, one could invoke the so called "Single 
state dominance hypothesis" 'M' according to 
which the sum in the eqs. ([6| and ([s]) is ex- 
hausted by its first term. The measured /3 
decay and EC ft values then make it possible 
to determine both the M'^" and M'^l;' . 
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Obviously, none of these methods is exact, but their 
combination has, perhaps, a chance of constraining 
the value of M'^^ substantially. Examples of appli- 
cation of the latter two items are shown in Table Ull 
That method can be used, obviously, only for the 
nuclei where the corresponding experimental data 
are available. 

Comparison of the NMEs Mj^j, and Mjl^ in Ta- 
ble |T] tells us, right away that, at least within the 
QRPA, the summation in the Eqs. Q and Q con- 
tains both positive and negative parts (see also Fig. 
[5]). This is obviously so since for most nuclei the 
quantity — {Mi -I- M/)/2 in Eq. ([s]) becomes 
negative, while each of the denominators in the Eq. 
([g]) is positive. Hence, we cannot expect good agree- 
ment between the M^^ from QRPA and those from 
the items 2. and 3. above. And, moreover, we 
cannot expect that SSD is a valid hypothesis for all 
candidate nuclei. Comparison of the corresponding 
entries in Tables [T| and [TI| confirms that expectation. 

Since there is a substantial experimental activity 
devoted to the determination of the strengths, it 
is worthwhile to examine in more detail the some- 
what unexpected finding that in many cases M^*^ 
and M'^i have opposite signs. Obviously, this has 
to do with the different weight of the corresponding 
terms in the Eqs. (|6| and We plot in Fig. [t] 
the corresponding running sums as a function of the 
excitation energy in the intermediate nucleus. One 
can see that the negative values of M'^^' arise from 
excitation energies Eg^ > 10 MeV that are difficult 
to explore experimentally. 

The negative contributions to M^" and M^l" from 
higher excitation energies cause in several nuclei 
even the reversal of the sign of M^l" to the nega- 
tive one. While, clearly, there is a substantial (3~ 
strength at these excitation energies, QRPA pre- 
dicts that there is a sufficient /3+ strength there as 
well, leading to the reduction of the M^'^ and M^l" 
visible in Fig. [7] Our QRPA calculations suggest 
that about 0.2 units of the B(GT) /3+ strength is 
distributed among states with Eg^ > 10 MeV in all 
considered nuclei. Such f3~^ strength has not been 
observed experimentally so far. It remains to be 
seen whether it exists at all, or is hidden in the 
"grass", i.e. distributed among many weak states 
that escape identification. Until this dilemma is re- 
solved we cannot decide whether the closure matrix 
elements M^l' in Table I are realistic or not. 

In the previous section we discussed the phe- 
nomenon of quenching of the axial current matrix 
elements. Fig. |6] suggests that using the effective 
g'^^ < 1.27 reduces the negative contribution of 
the higher lying 1+ states to the matrix element 
M'J'^. To see how large that effect might be we per- 
formed QRPA calculation with g'^-^ = 0.9 based on 
the empirical evidence that the degree of quenching 




FIG. 7: The running sums of Af^" (upper panel) and 
M^l' (lower panel) for selected nuclei, (color online). 
qa ~ 1.269 was used. 



increases with A. The resulting que nched matrix 



elements M'J'^ are shown in Table 



III 



While, 



marked earlier, it is unknown whether all mutipoles 
are affected by the axial current quenching, not only 
the GT 1+ states, we nevertheless show in the same 
Table the values of quenched Mq^ and of the full 
NME M'°'' with and without the closure approxi- 
mation. If quenching would not affect these matrix 
elements, their magnitude would be enhanced by the 
factor (1.269/0.9)^ ~ 2 making them substantially 
larger than the values in Table |l] 

Since our goal is the determination of the GT 
part of the Oi^PfS matrix element Mq^, a priori the 
knowledge of the M^i" , which depends only on the 
1+ intermediate states, is insufficient. According to 
the Eq. (11) we need for that purpose the func- 
tion C^;''(rythat depends, in principle, on all inter- 
mediate multipoles. However, is we could use the 
expansion of the spherical Bessel function jo{qr) in 
Eq.( 12 ) in powers of qr and keep just the first term, 
the neutrino potential H{r , E) would be represented 
by a constant and the Eq.( 11 ) would predict a sim- 
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TABLE II: The 2z^/3/3-decay closure nuclear matrix element \M'^i\ evaluated using the Single State Dominance 
Hypothesis (SSD) and with help of the measured strengths in charge exchange reactions (ChER). The adopted 
values of the 2i//3/3-decay half-times T'^J^"^^, taken from Ref. [17] are also shown. In the ChER case the matrix 

elements IM^"] and M^" have been determined by assuming equal phases for its each individual contribution. 



Nucleus [y] 


SSD 




ChER 




\M'^"\ [MeV-i] 




\M^"\ [MeV-i] 






4.4 X 10^^ 






0.083 


0.220 [25] 




1.5 X 10^^ 






0.159 


0.522 [26] 




2.3 X 10^^ 








0.222 [27] 




7.1 X 10^* 


0.208 


0.350 [29] 






"«Gd 


2.8 X 10^^ 


0.187 


0.349 [29] 


0.064 


0.305 [28] 




1.9 X 10^" 


0.019 


0.0327 29* 







TABLE III: The experimental values M^^p and the QRPA values of M'^i" evaluated with quenching at qa = 0.9. 
The corresponding quenched values of the Oi^/3/3 matrix elements are also shown. 

NME "'^Ge ^^ge ^'^Zr ^°°Mo "^Gd ^^^Te ""Te "^Xe 

2i/^^-decay NMEs 

|M^^;| [MeV-^] 0.136 0.095 0.090 0.231 0.126 0.046 0.033 (0,0.031) 
M^?" 0.336 0.100 -0.210 -0.205 0.179 -0.146 -0.169 (-0.25, -0.047) 

0i//3/3-decay quenched NMEs within closure approximation 
M'St-ci 2.50 2.10 1.20 2.23 1.63 2.09 1.77 (0.86, 1.08) 

M^?" 3.82 3.25 1.91 3.49 2.45 3.28 2.82 (1.43, 1.72) 

0!//3/3-decay quenched NMEs without closure approximation 
Mgt 2.59 2.20 1.33 2.45 1.71 2.23 1.91 (0.93, 1.14) 

M'°" 3.90 3.34 2.05 3.71 2.53 3.44 2.96 (1.51, 1.78) 



pie proportionality between M'^^ and M^'j.. How- 
ever, such an expansion does not work. In reality 
in Eq. ( 12 ) qr > 1 and we cannot approximate the 
neutrino potential H{r, E) by its value at r 0. 
Hence, we do not expect a proportionality between 
M'^i and Mg.^ and the QRPA evaluation supports 
this conclusion. 



responding nuclear matrix elements on the distance 
r between the pair of neutrons that is transformed 
in the /3/3 decay into a pair of protons. Thus (see 
Eqs. ((Sf, ^ and ^) 



Jo 



(16) 



VI. CONCLUSIONS 

Since the nuclear matrix elements M^'^ must be 
determined theoretically, it is of obvious interest to 
search for any relation between their numerical val- 
ues and other quantities that are either known from 
experiments or at least constrained by them. Here 
we describe such a relation between the dominant 
Gamow- Teller part M'^j. of M^" and the matrix el- 
ement M'^i of the observed 2i//3/3-decay evaluated 
in the closure approximation. The relation is based 
on the evaluation of the auxiliary functions CQ^(r) 
and C^l'{r) that describe the dependence of the cor- 



represents the required relation. 

However, while the matrix elements M^" and 
M^l' depend only on the transition strengths and 
energies of the 1+ virtual intermediate states (they 
are pure GT quantities), the function C^l'{r) gets 
contribution from all multipoles. Thus, the relation 
that we found is an indirect one; even if M^l'' would 
be precisely known, the evaluation of the function 
C^i'{r) requires additional nuclear theory input. 



Nevertheless, the relation in Eq. ( 11 ) allows us to 



obtain a better insight into the problem of the A and 
Z dependence of the matrix elements M^'^ and M"'^. 
While the known M^" have a strong shell depen- 
dence, the calculated M^" vary much less. Analysis 
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of the functions Cli{r) and C'Qrp[r) makes it pos- 
sible to better understand where this fundamental 
difference comes from. 

We show that, so far, the QRPA values of closure 
approximation M'^l' matrix elements do not agree 
well with the same quantities based on the measured 
/3~ and strength functions and on the assump- 
tion of coherence (i.e. same sign) of contributions of 
individual states. Until this discrepancy is resolved, 
it is difficult to employ M^f in order to constrain 
the magnitude of the matrix elements Mq'j.. 
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